NASA Technical Memorandum 104744 


Space Shuttle Entry Terminal Area 
Energy Management 


Thomas E. Moore 


November 1991 


(MAS A-TM- 104744) SPACE SHUTTLc ENTRY 
TERMINAL AREA ENERGY MANAGEMENT (NASA) ^ 

50 P CSCL 


G3/01 


N92-19930 


Line) as 
0057113 


NASA 




NASA Technical Memorandum 104744 


Space Shuttle Entry Terminal Area 
Energy Management 


Thomas E. Moore 

Lyndon B. Johnson Space Center 

Houston, Texas 


National Aeronautics and Space Administration 
Lyndon B. Johnson Space Center 
Houston, Texas 


November 1991 




CONTENTS 


Section 

1.0 

1.1 

1.2 

2.0 

2.1 

2.2 

2.3 

2.4 

2.5 

3.0 

3.1 

3.2 

3.3 

3.4 

3.5 

3.6 

3.7 

3.8 
3 : 

3.10 

3.11 

3.12 

3.13 

3.14 
3 5 

3.16 

3.17 

3.18 


Introduction 

Guidance Modifications 

Entry to Landing Guidance Functions 

Development History 

Energy to Attitude Control . 

Optional TAEM Targeting (GTT) 

Bailout Modes 

Smart Speedbrake and Bank Limit for 0122 

Theta Limits 

Detailed Description of Algorithms 

Technique of Constraints Limits 

Overview 

TAEM and GRTLS Executives - TGEXEC & GRE V EC .... 

Initialization - TQfNlT & GRINIT 

Heading Alignment* TGXHAC 

Navigation User Poityneter Processor - HAC Control 

Ground Track Predictor - GTP 

Computations - TGCOMP 

TAEM Transitions - TGTRAN 

TAEM Body Vertical Acceleration - TGNZC 

TAEM SpeedbraK© - TGSBC 

TAEM Bank - TGPHIC 

Overview of GRTLS Open Loop Guidance Phases 6 to 4 

GRTLS Transitions - GRTRN 

GRTLS Body Vertical Acceleration - GRNZC 

GRTLS Alpha Recovery - GRALPC 

GRTLS Speedbrake - GRSBC 

GRTLS Bank - GRPHIC 


Page 

.. 1 
.. 1 
.. 1 

.. 2 
.. 4 

.. 6 
.. 9 

.. 10 
.. 11 

.. 12 
.. 12 
.. 13 
.. 14 
... 15 
.. 15 
... 16 
17 

... 20 
23 

,.. 28 
... 36 
... 36 
... 41 
... 42 
... 44 
... 44 
... 45 
... 46 


ORIGINAL PAGE IS 
OF POOR QUALITY 


in 


PRECEDING PAGE BLANK NOT FILMED 



1.0 INTRODUCTION 

An historical account of the development for Shuttle’s Terminal Area Energy Management 
(TAEM), is presented. A derivation and explanation of logic and equations is provided as a 
supplement to the well documented guidance computation requirements contained within the 
official Functional Subsystem Software Requirements (FSSR) published by Rockwell for 
NASA. This FSSR contains the full set of equations and logic, whereas this document will 
address just certain areas for amplification. 


1.1 GUIDANCE MODIFICATIONS 

The abort Glide Return to Launch Site (GRTLS) is a high altitude extension of TAEM. 

TAEM was initially developed to do its primary function of energy management with an energy 
controller, but has since been modified to an altitude controller for better energy management 
In the presence of unknown winds. 

Optional TAEM Targeting (OTT) provides greater flexibility for pilot control of the ground 
track by providing the option for an overhead approach to the runway in addition to a straight- 
in approach. This modification has been in use since STS-5. 

A b ail out guidance mode was added to implement one of the recommendations of the Rogers 
Commission Report for crew capability to escape from an Orbiter in controlled subsonic 
gliding flight. The software changes to both TAEM and GRTLS, along with the associated 
escape pole, were the means of implementing this recommendation starting with STS-26. 

A pending modification, Change Request 89979, involves a smarter speedbrake for energy 
control and a smarter vehicle vertical acceleration (NZ) limiter that limits total NZ by 
limiting bank attitude. 

A potential modification, previously considered and still a potential change, is Theta Limits 
which would compute pitch attitude limits to constrain airspeed in the presence of degraded 
airdata. 


13 ENTRY TO LANDING GUIDANCE FUNCTIONS 

The Entry guidance flies the low lift-to-drag (L/D) and high angle-of-attack (a) lifting body 
vehicle from atmospheric entry to the higher L/D and low a aeroplane flight region of TAEM. 
The F ntry vertical component command of the aero lift vector commands the bank angle 
magn rhtJe to control the downrange component relative to the runway landing site. The 
later# Component command of the aero lift vector, commands the bank angle direction (roll 
revert#^) to control the crossrange component relative to the runway landing site, whenever 
the liberal deviation from the target exceeds a deadband level. As the TAEM target altitude is 
approached, entry guidance transitions from the high a of about 40° to the TAEM level of about 
10°. Speedbrake during entry is to a preprogrammed profile to assist attitude control and is not 
used for guidance. 

TAEM starts at a mach level of 2.5 at about 85,000 ft altitude. The lateral component command 
of the aero lift vector commands bank angle magnitude and sign to control crossrange relative 
to a prescribed or calculated (direct route) ground track toward the runway. The vertical 
channel basically controls energy and altitude with a. Speedbrake control during TAEM 
continues from entry to a preprogrammed profile until subsonic and then controls airspeed 
indirectly through dynamic pressure. 
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TAEM delivers the vehicle into the final approach plane of the runway and then delivers 
control to the approach and landing guidance at from 5000 to 10000 feet altitude. The 
approach guidance continues to fly the reference glide slope with a, the lateral track with bank, 
and the airspeed with speedbrakes. The landing guidance then flares from the steep reference 
to a shallow glide slope and then flares to the runway for touchdown. 


TAEM 



Entry to Landing Ground Track Sketch 

The GRTLS guidance is identical to TAEM, except that it starts it is TAEM phase at mach 3.2 
and has three guidance open loop phases with closed loop flight control of either angle-of- 
attack or acceleration in the region between external tank separation (about mach 6) and mach 
3.2. 



2.0 


DEVEI 




r HISTORY 


The three flight control variables given for the development of TAEM are speedbrake for drag 
modulation below mach .95, normal acceleration along the body axis Z (NZ) for longitudinal 
state control, commanded from the guidance to the flight control which controls a to get NZ, 
and vehicle bank angle or roll about the velocity vector for lateral state control or direction of 
flight. 


Flared Rudder 
Speedbrake 
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Flight Control Variables 


Energy is the longitudinal variable to be controlled. 

E = wh + mV 2 /2 

E/W = h + V 2 /(2g) = h + qbar/(pg) 

During the early 1970's, many organizations were Involved In the development of the TAEM 
guidance. This author. In JSC’s Avionics System Division, produced a Variable Entry Point 
guidance. Candidates from others were the Racetrack from JSC’s MPAD Division, the VORTAC 
from Charles Stark Draper Labs, the Cylinder from McDonnell Douglas, and the Spiral from 
Rockwell International. A study contract was given to MDAC/St. Louis to pick, or blend Into 
another form, the final guidance. The result was a “Hybrid" that Incorporated good features 
from all the candidates. 





Candidate Guidance Schemes 
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Simplicity was a major emphasis In the formulation of the Hybrid guidance. Instead of 
computing a ground track to satisfy the energy dissipation required from any initial condition 
to a runway landing, as was done with the VEP and Racetrack, the Hybrid would use a 
prescribed direct path to a cylinder fixed on the same side of the runway as the vehicle, nstead 
of using an offset target to provide an acceptable ground track in case of a high energy 
overshooting approach, as was done with the VORTAC and Cylinder, the Hybrid would just turn 
away (S-Tum) from the target cylinder (Heading Alignment Cylinder-HAC) for extreme high 
energy cases, until the energy entered an acceptable level to head for the HAC. Instead of a 
complex coordination of both vertical and lateral channels, as was done with the Spiral, the 
Hybrid would separate the channels: a for vertical, and bank for lateral control. 


2.1 ENERGY TO ATTITUDE CONTROL 

This baseline guidance, established in 1975, was then subjected to additional testing that found 
some poor performance and response characteristics. This baseline performed the major task 
for TAEM of energy control, by controlling the energy with speedbrakes, and by controlling the 
dynamic pressure with normal acceleration NZ. The first problem was that a firm control of 
energy is not desirable for a situation where the vehicle flies toward the HAC in a tailwind, and 
then turns 180" into a headwind. Energy control will bias altitude (h) below nominal and 
dynamic pressure (qbar) above nominal so that energy remains nominal. At the end of the 
turn, the low h and high qbar are the reverse of what Is required to fly into a headwind. 



Energy Loss Turning into a Headwind 


The second problem was the oscillatory response characteristics of controlling qbar with just 
NZ caused by qbar being a function also of drag in addition to lift (NZ). 
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NZ Control of Altitude or Qbar 

The vertical acceleration component shows close correlation with NZ, 

g»NZ=L/M Cos(a)+D/M Sin(a)-g Cos0 Cos<|> 
hdbldot = L/M Cos(<|)) Cos(t) - D/M Sin(?) - g 

i.e., exactly equal when a = -yand <|> = 0, and therefore good predictable quadratic response 
can be expected from: 

NZC = ki (he - h) - k 2 hdot 

but, looking at qbar, 

qbar = p V 2 /2 

where, 

p = PO e‘ Kh 
pdot/p = -K hdot 
Vdot = -g Sin(y) -D/M 

and the derivative of qbar, 

qbardot = 2 qbar Vdot / V + qbar pdot/p 
= -2 qbar ( g Sin(>) + D/M) / V - qbar K hdot 

shows a gooc' correlation with hdot, which by itself, would give a good correlation to NZ by the 
above argument. But it also has a predominant drag term, where the representative numeric 
example: 

g=-12, g=32, D/M=2*32/4=16, K=5E-5, 
qbar=200, hdot=-150 

qbardot = +3 -8 +1 .5 

shows the -8 drag term to be larger and opposite in direction, and therefore poor unpredictable 
response can be expect from: 

NZC = -k 3 (qbarc - qbar) + k 4 qbardot 
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The change to controlling just the altitude component of energy with NZ solved both of these 
problems. Maintaining reference h eliminates erroneous altitude biasing of energy control, 
and the close relationship of NZ to altitude acceleration produces predictable and smooth 
response. For simplicity sake, the speedbrake was changed to the other component of energy 
which is velocity or qbar t and for the sake of approximately matching the controller of the 
Approach & Landing guidance, which is based on airspeed. 

This altitude controller also maintains control over both energy and dynamic pressure by the 
utilization of midvalue limits that are calculated. This will be discussed In more detail later. 


2J2 OPTIONAL TAEM TARGETING (OTT) 

The Shuttle flight crew in 1976 expressed an interest in an overhead (OH) type approach as an 
option to the straight-in (SI) approach to the HAC. The objective is to conserve energy for 
subsonic dissipation in the vicinity of the runway. Another reason which became the real 
selling point for making the change was thunderstorm avoidance at the Cape. The altitude 
difference between SI and OH increases the probability for successfully flying either under or 
over a weather disturbance along the landing site approach. 


Weather 

Disturbance 


Note: OH/SI toggling should be done eariy in 
Entry guidance to minimize trajectory maneuver. 



Range-to-go OH is 
14n.mi.>than SI, 

& altitude is about 
25000 ft higher. 


OTT Option 


The development of this guidance modification started In 1976 as the Reverse Targeting 
Scheme and later became known as Optional TAEM Targeting (OTT). The initial developers 
and testers were the author of this paper, Vance Brand/Astronaut, Bob McNenny/MDAC, Gil 
Carman/ JSC-MPAD, Chuck Bowser/JSC-Crew Systems, and Ellis Henry /JSC-MPAD. The 
Initial development was on an off-line Univac batch computer, and then evolved to more 
detailed testing and development on the Crew Procedures Simulator-Bldg 5, Shuttle Procedures 
Simulator- Bldg 35, Shuttle Engineering Simulator-Bldg 16, and the Shuttle Avionics 
Integration Lab-Bldg 16. The first flight of OTT was on STS-5 In 1982. 

In addition to the option of SI or OH, OTT added or modified the following features: 1. The 
Heading Alignment Circle-HAC was changed to a spiral. In three dimensions the circle is a 
cylinder, and the spiral a cone, and therefore the HAC is now a Heading Alignment Cone. 
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Spiral HAC 

RTURN = RF + R-|'P+ R2'P 2 
Range = (RF + Ri^/2 + R 2 'F 2 /3) v F/57. 

RF =5K to 14K, R1 =0., R2=.093 

2, The HAC radius is now adjustable after the \C turn phase starts. A low energy condition 
of sufficient magnitude will start shrinking the radius. 



Spiral Adjustment 


RF = RF -0.8 ( HREFOH - H) / ( 3H/3R • 3R/3RF) 
where, 

3H/3R = Tan(15°) 

3R/3RF = 'FV 57.3 
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3. Provisions for HAC turns greater than 360° were made because Entry guidance roll reversals 
can wrap an OH which starts less than 360° to greater than 360°. 



4. An energy dump phase on approach to the HAC was added to alleviate the ground track overshoot that 
results from large turn, supersonic HAC starts. Supersonic bank is limited to keep sonic boom overpressure 
low. This procedure biases the target energy lower then nominal for large turn angles, which produces a dive 
to dump energy, and then just prior to transonic, a pull-up maneuver targets for subsonic qbarmin. This 
converts kinetic to potential energy, allowing the HAC to be followed more closely; because the subsonic 
region starts sooner, the higher subsonic bank is used, and the velocity to be turned is smaller. 



\ 


Dive to Energy 
Bias 


Large Turn, High Energy HAC Overshoot 
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Energy Dump for High Speed Turn 

5. Also this low supersonic bank limit could, with default airdata, produce a serious HAC 
overshoot and loss of vehicle because of delayed transonic indication from AD. The solution 
implemented is that the low bank limit is increased if airdata is default. 


2.3 BAILOUT MODES 

The bailout mode was developed to provide an automatic flight system to maintain the Orbiter 
1 controlled flight to allow time for the crew to escape using the newly developed parachute/ 

, ole escape system. The design requirements: 

Airloads minimized by minimizing airspeed. 

Two fault tolerance on mode engagement resulting from either crew action or hardware failure. 

Nominal manual flight operations maintained. 

Angle-of-attack constrained to airdata maximum of 20*. 

Minimize software to fit within the API 01 B flight computer. 

The level B changes were to set a bailout flag if mach <.95 and the pilot had both moved the 
abort selection switch to the Abort to Orbit position and pushed the abort PBI. The bailout flag 
would then force GKTLS to re-ini tialize to the acquisition phase. 

The TAEM/GRTLS guidances monitor the bailout flag and the manual/auto status of both pitch 
and roll, so that command values are obtained at the point that the pilot has selected those 
values by going from manual to auto, individually in each axis. The pitch command selected is 
dynamic pressure, which is controlled with the standard logic, by making the minimum and 
maximum values the same as the command value. Also a new angle-of-attack pitch axis 
constraint is added for both TAEM and GRTLS. The roll command is simply the bank angle 
selected at the above snapshot. The speedbrake command is simply zero when the bailout flag 
is on. 
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2.4 SMART SPEEDBRAKE AND BANK LIMIT FOR 0122 


The baseline speedbrake controller has gone from energy to qbar and now eventually back to 
energy although with a different formulation from the original. The change request CR89979 
is scheduled for 0122, Jan 1993. 



* Equivalent signal : -EDTERI - (dEc / dR - dE / dR) dR / dt 


** Eilat 


tQ — { 


EDSMGAN 


Hr 


Smart Speedbrake 


The present qbar algorithm generally moves the brakes in the proper direction, except in 
transient conditions such as a vehicle (energy & altitude) high, qbar low case. Energy will put 
them out, whereas 0bar will tuck them erroneously in, until the altitude channel drives the 
speed up. 

The NZ and bank limit coordination are also Included in CR89979. This will produce a 
cooperative coupling, such as for a vehicle high condition during a turn maneuver in a 
crosswind. 
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Equilibrium component of normal acceleration, 
NZC Cos $ = Gravity (1 .0) Cos 0 


Total normal acceleration, 

Cfcslrod limit ± 2.2 [ NZ tota | ] _ CoS 0/ CoS<|) + [ NZC ] M± - 5 


To constrain NZ total to NZmax 
then solve for bank limit, 
NZmax - Cos 0 / Cost)) lim + NZC 
Cost)) lim - Cos 0/ (NZmax - NZC) 
t)i lim « Cos-1 (Cos 0 lim) 


Smart NZ - Bank Limiter 


The new computed baf)|( |jmit will provide more lateral acceleration to counter the crosswind 
as the vehicle pitches down toward the reference altitude. The NZ Is 1 . ited for structural 
reasons, and now If the vertical channel does not require much lift, then the lateral channel 
can get more force up to the NZ limit. The reverse situation can also be found where the vertical 
channel can get the maximum from NZ if the lateral does not require much force. 


2.5 THETA LIMITS 

Erroneous airdata to TAEM could result in the NZ(qbar) control laws either overspeed 
overloading the vehicle, or underspeed stalling, or losing energy on the backside of the L/D 
curve. This can be prevented with additional constraints on NZC of upper and lower pitch 
limits (0). 


6max = f(V, qbarmin, Speedbrake, Bank) 

0min = f(V, qbarmax, Speedbrake, Bank) 

These limits can be determined empirically by flying an open loop entry simulation to the min 
or max qbar for various setting on speedbrake and bank. This has been done and the results 
were applied to the pilot displays to allow a manual override of the guidance. This could be 
applied to the guidance utilizing another midval limiter added to the present NZC. 
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Potential Theta (6) Limits for TAEM 

An open Issue though is where to place this logic. It is a guidance function, but it may require 
the higher computation frequency of the flight control loops. 


aO DETAILED DESCRIPTION OF ALGORITHMS 


ai TECHNIQUE OF CONSTRAINTS LIMITS 

A utility routine MIDVAL is used throughout the guidance that selects the middle value of three 
inputs that can be mixtures of variables and constants or all variables. A cascade of these 
software devices has the effect of prioritizing constraints. For example 



The lowest level constraints say that 12 will never exceed II , nor go below 13 , assuming II > 13. 
The highest priority constraints say that no matter what 1 1 ,2 or 3 do, 0 will never go above K1 , 
nor below K2, assuming K1 > K2 . 

The most detailed application of this technique is seen in the TGNZC routine of TAEM where 
three midvals are cascaded. 
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TGNZC Constraint Priorities 


with the result: 

1 . NZ(h) is the nominal driver, so that usually NZC = NZ(H) 

2. But if energy goes too far beyond energy nominal then the dynamic commands 
NZ(E/Wnom+8000) or NZ(E/Wnom-4000) take over. For example, a simple straight-in approach 
to the runway in a headwind, could possibly not make H to the runway if the vehicle starts pitching 
down to a less efficient L/D as soon as h is satisfied. 

At V = 2000, where 


E/W = h + V * V / (2g) = h + qbar / (p* g ) 

AE/W = Ah + Aqbar / ( p * g) « Ah + V*V * Aqbar / (2g*qbar) 

AE/W = Ah + 500 Aqbar 

Aqbar can be about 80 low(140) from the nominal(220) for a L/Dmax case, which then at mach 2 
would be equivalent to 40,000 ft of altitude. That is E/Wnom is not achieved until h goes 40000 
higher than its nominal. 

3. But regardless of the h or E/W situation, a higher priority is NZ(qbar), i.e., 140 min for max L/D or 
300 max for max dive. If energy were low, it would only go lower if qbar were allowed to go below 
qbarmin. 

4. And the highest priority is vertical acceleration which prevents the wings from being pulled off 
trying to get to an optimum qbar. 


3J2 OVERVIEW 

The remaining discussion -ill follow the format of the FSSR’s for both TAEM and GRTLS. 

There are unique routine r or GRTLS which are given names that start with GR . The 

basically common routine lsed by both TAEM and ORTLS are given names that start with 
TG The basic functions of TG routines between TAEM and GRTLS are the same, although 
they may not be exactly Identical. For example. TGCOMP has more linear reference profile 
segments for GRTLS than TAEM. 

Any I-Load constants values should be considered approximate where exact values should be 
obtained from the I-Load database. 
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The logic flow: 


CONDITION \ 

then j 

Block 1 

1 /I 

else 


i 

V 


Block 2 


means IF (CONDITION) is true, then do Block 1 , 
else, if false, do Block 2. Return and continue past 
CONDITION when finished either. 


3.3 TAEM AND GRTLS EXECUTIVES - TGEXEC & GREXEC 

The flow chart of all TAEM subroutines also shows the Interface with the pilot for OTT, the 
various phases of TAEM, and a brief description of changes to SB & Bank with CR89979. 



TAEM Executive Row Logic - TGEXEC 

The flow chart of all GRTLS subroutines also shows the Interface with the pilot for OTT, and 
the various phases of GRTLS. 
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GRTLS Executive Flow Logic - GREXEC 


3.4 INITIALIZATION - TGINIT & GRINIT 

All variables t v ;>■ require inlt -alization are done here, for either the first pass through TAEM 
or GRTLS, or wnenever during these guidances, the opposite end of the runway or a new runway 
is selected. 


3J5 HEADING ALIGNMENT- TGXHAC 

The inplane X coordinate points of the Heading Alignment Cone are computed from the 
corresponding pre-flight trajectory design input (I-Load) altitude points. 
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Geometry calculations of the Heading Alignment Cone 
-TGXHAC 


The vehicle nominally arrives In the final approach plane on the steep glide slope at HFTC 
prior to the Approach & Landing Guidance at HALI, but If an energy low condition were 
encountered and the pilot assesses that HALI will not be achieved, then the pilot can move the 
HAC to XMEP to shorten the range to the runway. 

The adjustment of the glideslope intersect point XA (IGI) is an autoland adjustment for 
headwind. The adjustment should be made in TAEM, to allow time to settle on the new 
trajectory, but it is not intended to be a TAEM energy control variable. 


3.6 NAVIGATION USER PARAMETER PROCESSOR - HAC CONTROL 

The pilot can toggle the HAC between overhead (OH) and straight-in (SI) with a keyboard item 
entry. The software is outside the guidance and documented in the Navigation FSSR 6/30/85 
pages 4-207 to 210. 



NAV UPP HAC Control 
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The HAC Information sent to the guidance specifies the side of the runway for the HA and the 
targeting for either OH or SI. The HAC turn angle PSHA is reset to PSHARS at a toggle to 
restart properly from a potential wrap around condition where PSHA > 3C 


&7 GROUND TRACK PREDICTOR - GTP 

The ground track, or trajectory component in the ground plane, is predicted for three major 
phases by calculating the turn geometry of the acquisition phase by rotating the present 
velocity vector toward the final turn spiral HAC, and then wings level flight to the HAC. The 
HAC phase distance is calculated from the acquisition point around the HAC and into the final 
approach plane, and then to the runway for the pre-final phase. 



The GTP is not concerned with energy state at the present or end conditions. Other elements, to 
be discussed later, will guide the vehicle energy to a reference at the predicted range. 

The navigation state of the vehicle and the state of the HAC are specified relative to the runway 
from which the GTP computes the state of the vehicle relative to the HAC and the vehicle range 
from the runway. 

Working backwards along a trajectory, the last part of GTP is simply the direct range vector to 
the runway. 



Prefinal Approach Range Prediction 


17 


This direct range algorithm used to start with the last guidance phase (3), but as part of the OTT 
modifications, this start was delayed until XCIR < DR4(2000) to provide a smooth range 
calculation transition. A smooth bank calculation transition occurs with the phase 3 point at 
DR3(5300 ft). 

The next part of GTP Is the analytical range around the HAC. For the spiral HAC, 

RTURN = RF + Ri PSHA(deg) + R2 PSHA 2 
(14K) (0) (.093) 

The range around the spiral, 

RPRED2 = J Tan Vel • dt = / (RTURN*PSHAdot/57.3) *dt 

PSHA 

= jRTURN/57.3 *dPSH A 
0 

= (RF-PSHA + R r PSHA 2 /2 + R 2 *PSHA3/3)/57.3 


and to the runway. 


RPRED2 = RPRED2 + XHAC 



The acquisition phase establishes RTAN the same as the HAC phase. The velocity vector is 
then projected through a circular tum until parallel to RTAN. From that point the RC vector is 
projected to the HAC Intersect. 
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ARCAC - RTAC • | DPSAC | 

RPRED - ARCAC + RC + RPRED2 

Acquisition Range Prediction 

The circular turn radius is calculated from the present velocity and an average bank of mach. 



30 . 




Mach 
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For ydot =0, 


L/M Cos(<(>) = g Cos(t) 


RTAC = VH 2 /(L7M Sin(<{»)) = VH 2 /(g Cos(y) Tan(<J>)) 
= VH»V/(g Tan(<)>)) 

Acquisition Turn 


An additional terni is added for GRTLS because of the early S-Tum at high velocity with a low 
bank rate of only 5 deg/ sec. 


v 



Acquisition Range Prediction 
Additional Term for GRTLS 


&8 COMPUTATIONS - TGCOMP 

General calculations that will be used elsewhere are basically done here, such as energy, 
reference values for energy, altitude and dynamic pressure, and boundary values for energy. A 
few exceptions, though, are that the action of filtering qbar and the control of the HAC final 
turn radius are done here. 

There are three variable indexing parameters for TAEM.: 1. IGI, which was discussed in 

TGHAC, has two settings for glide slope intersect for final approach; 2. IGS has two settings to 
allow reshaping of profiles for an optional heavy payload return from orbit. GRTLS does not 
have this index, because the weight is known before lift-off; 3. IEL allows two segments for 
TAEM and four segments for GRTLS to generate reference profiles. 
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Energy Reference (Nominal) Profile 
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GRTLS also does additional indexing with IES, IEST and IEM to generate S-Turn and 
minimum energy profiles. These indexing differences are the main reason that the TGCOMP 
software, used by both TAEM and GRTLS, are not identical, although the functions are the 
same. 

The objective of the Energy Dump Maneuver, in cooperation with the Pull-up Maneuver in the 
TGNZC, is to target for subsonic conditions at HAC initiation, to facilitate tracking the HAC. 



Energy Dump 

The empirically determined R2MAX is the maximum value of RPRED2 to start the turn phase 
at the nominal energy and be at a low enough speed to track the HAC. The empirically 
determined ESHFMX is set to minimize energy loss for large turns > 400". 

The energy profiles are computed here and used for decisions in TGTRAN, or constraints in 
TGNZC. One exception is that EMOH for TAEM, but not GRTLS, is computed in TGTRAN. 



Energy Profiles 
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The energy boundary lines of ES and EMEP represent the limits of wings level dive or range 
stretch respectively. The max stretch for an overhead approach, EMOH, is less than EMEP 
because higher energy dissipation during the HAC turn. The EMAX & EMIN are used to 
constrain the altitude controller In TGNZC. Energy Is a better control parameter for a straight 
in approach with either a wind condition where energy will bias altitude to compensate 
velocity, or a navigation altitude error. But altitude is a better control parameter for non- 
straight in wind conditions. The energy constraint on altitude tends to capture the advantages 
of each. The bias values from EN are empirically determined from 180" turns into either a 
head or tail wind on final approach. 

The altitude reference is computed from linear and cubic function segments. 



Altitude Reference Profile 

The equations for the Tan of the flight path angle reference, DHDRRF, is obtained by 
differentiating the href equations. 

The dynamic pressure reference is generated from two limited linear segments. 


qbar 



Dynamic Pressure Reference Profile 


The radius of the HAC spiral at the final approach, RF, is adjusted if the altitude goes below 
HREFOH during the turn phase if the turn angle > 90". 


Altitude 


HREF 

HREFOH 

h 



Range 


Spiral Adjustment 
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The adjustment utilizes a filter. 



RF Filter 

Assuming the HAC flight path angle Is about 15', then 

h = RTan15 , he = RcTan15 
R = RF • ¥ , Rc = RFC • W 

Solving for command values, 

RF c - RF = (h-hc) / (¥ Tan15) 

and substituting HREFOH for command altitude for the rate of change of RF, 

DRF=(-.8/Tan 1 5) (H R EFOH-h)/(PSHA/57.3) 

The Input dynamic pres' ire to the qbar filter. 



Qbar Filter 


Is selected from either measured airdata. If it Is good and V < 1500, or from navigation derived 
qbar. 


09 TAEM TRANSITIONS - TGTRAN 

This routine, used by both TAEM and GRTLS, performs transition of guidance phases and 
issuance of alert conditio >«. There are additional GRTLS phase transitions In GRTRN. 
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IPHASE 

Description 

0 

S-Tum 

1 

HAC Acquisition 

2 

HAC Turn 

3 

Pre-final Approach 

4 

Alpha Transition 

5 

NZ Hold 

6 

Alpha Recovery 

End 

Start Approach & 

TAEM 

Landing Guidance 


TGTRAN 


GRTRN 


In addition, the bailout mode is managed. When the bailout flag has been set then individually 
the pitch and bank flags are set for bailout control of each axis whenever auto mode of that axis 
has been selected after having been in manual mode. The manual mode drives to the 
conditions that will be held by the auto mode. 

The transition out of TAEM and into AL guidance for normal non-bailout operation is a 
function of four error functions on h, Y, y and qbar, or it is forced at h of 5000 ft. The 
transition criteria was provided for TAEM implementation by the AL development group. 

The transition to phase 3 usually occurs at a range of RPRED3, but will be forced to phase 3 at h 
< HMIN3(7000) to enable the AL transition for a low energy approach which can occur only 
from phase 3. 

The transition to phase 2 occurs when the vehicle gets to within 10% of the turn radius. 

RCIR < 1.1 RTURN 



This empirical transition point provides good bank transition by avoiding the HAC overshoot 
of a later transition and a bank command sign reversal of an earlier transition. This earlier 
start problem results from the radial rate damping term for phase 2 bank command. The ideal 
RCIR start produces and Initial bank command equal to the steady state value required for the 
turn. 
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Initial HAC Bank Command 


The transition to an S-Tum is done if E/W exceeds the profile line ES. The transition from the 
S-Tum to phase 1 occurs in TAEM when the energy returns below ES by a bias amount ENBIAS, 
and in GRTLS when energy goes below the profile EST(see TGCOMP). The objective of the S- 
Tum is to do an open loop maneuver turning away from the straight- in solution until energy 
gets to an acceptable value(EST) for straight-in. Attention is also given to the direction of turn 
because a turn which intersects the HAC in the wrong direction is undesirable. 



The direction of turn logic that unwraps the HAC during an S-Tum prevents this problem, 

S (Sign of Bank, + right) = - YSGN 



t 
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Right Direction S-Tums 

Additional turn direction considerations are given for turns < 90'. 



To avoid an S-Tum geometry problem for low range condition where an acquisition solution is 
not achieved after the S-Tum, 
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S-Tums are inhibited if DRPRED < RMINST (20n.mi.). At this minimum range S-Tum, 



there is sufficient room to turn and acquire the HAC. 

The S-Tum is also Inhibited in TAEM for turn angles > PSSTRN(200'), 



because the high energy HAC overshoot will gain enough range and transition to a satisfactory 
energy-range state without encountering a geometry problem. 

But for GRTLS that nominally requ iVgs S-Tums in phase 4, PSSTRN is set at 1000' to enable S- 
Tums. 


Low energy alerts are issued to tjje Allot for his action. An energy lower than EMOH during 
overhead approach is suggesting -fpflepilot that he consider downmoding to straight in. 


ORIGINAL PAGE 
OF POOR QUALI1 
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An energy lower than EMEP is advising that the HAC be moved closer to the runway. 

aiO TAEM BODY VERTICAL ACCELERATION - TGNZC 

The normal acceleration command from guidance to the flight control is an incremental 
command of acceleration along the body negative Z axis. 



and the flight control system adds to this a feed forward, steady state gravity compensation 
term computed by matching the vertical component of NZ with gravity along that direction. 



iiiiiiiiHiiii 111 uni 111 1 



Flight Control Nz Loop 


The primary guidance vertical controller on altitude. 



14K 24K 


Altitude Controller 
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is simply a second order controller, with the assumption that DNZC produces hdbldot without 
time delay or error. The stability and response characteristics are specified by design to be at 

a natural frequency con, and a damping ratio £, determined by the products of gains and the 

ratios of gains. 

co n = GDHV(HDREQG *.322) 

= .054 to .179 

The gain ratio is held constant as frequency varies with altitude and therefore the damping 
ratio is slightly under critical damping. 

4C 2 = .322 GDH / (HDREQG*GDH) = 3.22 
C=-9 


A more rigorous evaluation of the altitude controller transfer function involves NZc to hdbldot 
transfer. 

NZc = hdbldotc / g 

Assuming that Az converges to the command state, i.e., flight control gains » guidance gains 
then 

Az = NZc + Cos 9/ Cos c|> 

The measured acceleration is given by, 

Az = L (Cosa+ Sina/(LVD))/(Mg) 
and altitude acceleration is given by, 

hdbldot = L (Cos <)> Cosy - Siny(L/D))/M - g 

Solving from the above, 

hdbldot= Q ( hdbldotc + g Cos 0 / Cos 4> ) - g 

where, 

Q = ( Cos $ Cos y - Sin y / (UD)) / ( Cos a + Sin a /(L7D)) 

Substituting hdbldotc = ( CDn/(2Q ( he - h ) - hdot ) 2£(Dn , and transforming to LaPlace, with the 
linearizing assumption that Q is constant, 

h(S)=(Q CD n 2 hc(S)+L{g(QCos 9/Cos <(> -1)}) 

/ ( S 2 + Q2C(0 n S + Q (On 2 ) 

Numerical evaluation with L/D=4, <x=10, y^-10, 9=0, and <j>=0 gives, 

Q= 1 

h converges to he without error 
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With <t> =45, and COn =.179, 


Q = .72 

h converges to he +he 
he = g ( Q Cos9/Cos<|> - 1)/(Q © n 2 ) = 26ft 
and the actual damping decreases 

C actual = 'Tq •£ design (.9)= .76 

A modification that divided hdbldotc by Cos <f» would change this damping decrease to a small increase, 
.92. 

The guidance has been formulated without this Cos<t> division and, so far, no justification has been found 
for a change. 

Stability analyses of Gain and Pbnse Margin, Root Locus, and mainly integrated simulation 
were conducted for the total G&C ’em where flight control gain is not » guidance gain. The 
next lowest gain(qc/Aa) in the fb fc control NZ loop is .47 at qbar of 200 which is not much 
higher than the highest guidance gain of .322. The stability analyses, though, have all verified 
satisfactory response characteristics and control margin. 

The secondary vertical controller' ire for the dynamic constraints of dynamic pressure and 
energy with NZ. 

The maximum dynamic pressure constraint QBMXNZ is determined from hinge moment 
constraints. The minimum QBMNNZ is the results from a flight at maximum L/D and is a 
function of weight and bank. The profiles are generated as a function of mach. 


340 

QBMXNZ 


275 ( 300 for GRTLS) 



Dynamic Pressure Constraints 


To facilitate trajectory tracking of ie HAC, a pull-up maneuver that lowers QBMXNZ acts in 
cooperation with the energy dump of TGCOMP to target for subsonic initial conditions for a 
large turn HAC. The maximum turn angle at which the HAC initiation is subsonic, with a 
nominal energy approach, is approximately 270" at a range-to-go RPRED2 of R2MAX. At 
higher turn angles the supersonic low L/D and low sonic boom bank limit can cause the range 
to increase at HAC initiation as the vehicle flies past the HAC. 
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High Speed Turn 


The energy dump shifts the reference energy profile to target the same energy state of R2MAX at 
RPRED2. The pull-up starts at an energy state EQLOWU 85,000ft, by lowering QBMXNZ by an 
amount depending on the range error term and any flight deviation from the shifted energy 
reference. 



Management of HAC Initiation Energy State 

This energy dump & pull-up maneuver lowers the energy deficit that must be made up 
subsonically. 

The secondary guidance vertical controllers for minimum and maximum dynamic pressure 
constraints each use the same flow diagram. 
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QBMNNZ 

& 

QBMXNZ 


.075 .0625 


- QBG1 

1— 

-GBG2 - 






QBNZUL 

& 

QBNZLL 


QBARF 


QBD 


Dynamic Pressure Controller 


An added feature In GRTLS Is that QBG2 Is a function of mach. 

Analytical stabili analysis of qbar is not as clean as with altitude because there is no simple 
relationship with NZ, but is a function of many variables and derivatives. Also, it is non- 
linear. 


qbardbldot=pVVdbldot+pVdot 2 

+2pdotVVdot+pdbldotV 2 /2 

The alternative approach used was to determine gains and stability empirically. The response 
at the flight boundaries of either qbar min or max, is satisfactory. The unpredictable and 
oscillatory response, though when used as a primary controller, was a part of the reason for 
selecting altitude as previously discussed under “Energy to Altitude Control." 

The secondary guidance vertical controllers for upper and lower energy constraints each use 
basically the same controller as altitude, except that energy error replaces altitude error, and 
altitude rate dai ng remains the same. The gains were given new names, but to date, the 
values have beer. same. 


VH 


EMAX 

or 

EMIN 


iDHDRRF | 


iHDREF 


.1 


-o 


QEUL 
or ■ GDH 
GELL 


E/W 


QH J22-G0H 


HDOT 


DNZC 


1/G 


GEHDUL 


or 

GEHDLL 
-,322/g- .01 

Energy Controller 


This controller assigns all of the energy change to the potential component. The final energy 
will converge to the appropriate potential (initially unknown because the Integral of drag*time 
is the only thing that changes energy) and kinetic components. The stability is approximately 
that of the altitude controller, but is mainly verified by simulation to be satisfactory. 

Constraints on angle-of-attack are imposed for GRTLS only. 
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20 



.5 1 . 1.5 2 . 2.5 3 . 35 

Mach 


Alpha Constraints 


The controller for a constraints, 


ALPUL 

OR 

ALPLL 



Alpha Controller 


The convergence to a steady state will occur when ydot = 0. 

ydot = (L/M Cos<|> -g Cost) / V = 0. 
NZ=(Cosa+Sina/(L/D))L/{Mg)=DNZC+Cos0/Cos<]> 

Solving for DNZC and ot£ = ac - a, 

DNZC=(Ca+Sa/(L/D))Cy/C<|> - C0/C0 = ae 25.*.005 
Numerical evaluation with 

L/D=4, a=10, t=- 10, 0=0, and 0=0 , and 
1 5 20 0 20 0, gives, 

as = .1 to 1.8* 


and therefore a converges closely to ac. 

The baseline altitude and the six other dynamic constraint signals that have been generated 
above are now combined using the "Technique of Constraints Limits” that was previously 
discussed. 
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Vertical Controller Constraints and Filter 


where the alpha constraints are for GRTLS only. 

The bailout mode takes a snapshot of dynamic pressure QBARS whenever bailout_pltch goes 
from false to true. The upper and lower qbar constraints axe then set to QBARS so that 
NZ(QBARS) will be midvalue selected, regardless of NZ(h, E/W or alpha), to go to the CQG filter. 
An additional angle-of-attack constraint is added for both TAEM and GRTLS (which replaces 
ALPNZUL previously computed in GRTLS). 



Vertical Controller Constraints for Bailout 

In the pre-final approach phase 3(and not bailout), when it is too late to do anything about 
energy and maximum control tightness is desired, the filter is removed and Just qbar 
constraints are used. 


QBNZUU 

M 

D 

DNZUL 

M 

1 

D 

NZC 

DNZC 

NZC 


V 


V 

to 


A 


A 

FCS 

QBNZLL 

L 

DNZll 

L 



Vertical Controller Constraints for pre-final 
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3.11 TAEM SPEEDBRAKE • TGSBC 




The speedbrake controller is basically a proportional plus Integral on qbar error. An energy 
term has been added to compensate for situations where qbar and energy may not be In sync, 
such as, a 180" turn from tall to head wind. The energy can drop sooner than qbar. The energy 
term will start to modify the brake upper limit at -427 ft., and will have It reduced to full in by - 
7000ft of energy error. A replacement of this controller with a direct function of energy, at 
least for the early part of TAEM, was discussed previously under “Smart Speedbrake for 0122." 

The bailout mode sets the speedbrake command DSBC to zero so that the output DSBCAT will 
be the lower limit DSBCLL function of mach. 


3.12 TAEM BANK - TGPHIC 

The upper and lower limit Imposed on all the bank commands at the end of this routine are 
computed at the beginning as functions of mach and phase. 


PHILIMIT 



Bank Limits 
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The S-Tum bank command is simply at the bank limit, in the direction of "S", as determined 
in TGTRAN. 

PHIC = S • PHILIMIT 



The acquisition bank command is proportional to heading error. 

PHIC = GPHI (2.5) • DPSAC 



Acquisition Phase 1 Bank Command 
or Phase 2 with RERRC > RERRLM (7000) 


The relationship of lateral acceleration and bank with steady state vertical conditions, 

L/M Cos <(* = g Cos y 
L/M Sin <{> = a 
Tan<|> = a/ (g Cosy) 

<j>’ ~ 57.3 Tan <|> 

The rotation rate of the velocity in degrees per second, 

© = a N = g Cos yTan 4> 57.3 / V 
4>* ~ to V 57.3 / ( g Cos Y57.3) 


Substituting command values. 


<J>*c ~ toe V / (g Cos y) 
= 2.5 DPSAC 
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Solving for the control gain, 


coc / DPSAC = 2.5 gCosy/V 

This lateral controller gain Is small compared to the next Inner loop gain(.5 to 1.2) In the flight 
control, 



Flight Control 

Lateral Acquisition Controller 


and therefore the response Is exponential at about 1/.08 = 12 second time constant. 

The HAC turn bank command Is proportional to position and rate errors relative to the HAC, 
and includes a feed forward centrifugal force acceleration term. 

PHIC = -YSGN(RdbldOtref + GR*AR+ GRDOT »ARdot) 



The relationship of lateral acceleration and bank with steady state vertical conditions. 

Tan 4>c - Rdbldotc / (g Cos y) 

and the partial derivative, 

3<J>*/dRdbldot = 57.3 Cos 2 <|> / (g Cosy) 
are used to derive the lateral controller. 
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Rdotrf 


RTURN 


RCIR 


Rdbldotrf « 


GR/GRDOT 

h? 

■ GRDOT/K 



.006/2 ^ 
.025 
(fln/2£ 


57.3 Cos 457(g Cos H -YSGN 


st/3 

.22 

2CC0n 


K 

.9 


Lateral HAC Turn Controller 


which has response characteristics for a natural frequency of .074 f and is overdamped with 
damping ratio £of 1.5. 

The reference radial rate is derived by differentiating the spiral radial command. 


RTURN = RF + R1-PSHA + R2*PSHA 2 
RDOT = (R1 + 2 R2*PSHA)*PSHAdot 

PSHAdot = -VH 57.3 / RTURN 
RDOTRF = -VH R1+ 2 R2*PSHA) 57.3 / RTURN 


The reference radial acceleration counter balances centrifugal force, 

Rdbldotf = -Vt 2 / RTURN 

The Implementation of this term involves the conversion to <|>c, 

PHIP2C = Rdbldotf 57.3 / g 
and then a quadratic > urve fit from Tan4>c to <|>c. 

Tan<f>c = PHIP2C 

PHIP2C= PHI2C1*Tan«|)C - PHI2C2*Tan 2 <(>c 
1.13 -.0055 

4>C PHIP2C 

30 31.4 

45 46.7 

60 58.0 

The pre-final bank command is proportional to lateral position and rate errors relative to the 
runway. 


PHIC = - GY-Y - GYDOT -Ydot 


Y 


Pre-final Phase 3 Bank Command 
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The relationship of lateral acceleration and bank with steady state vertical conditions, 


Tan <j>c = Ydbldotc / (g Cos y) 
<)>c ~ 57.3 Ydbldotc / (g Cosy) 


The lateral controller. 



(On/2; 2;®n 


Lateral Pre-final Controller 


has response characteristics for a natural frequency of . 15, and a damping ratio £ of .73. 

A fader which starts at the first pre-final pass with ISR=5 allows the bank command to change 
by only 1/5 the difference between old and new command. On the second pass 1/4 is dumped in, 
and so on, until at ISR=1 the system fully uses the phase 3 PHIC. 

Simulation testing revealed two oscillation problems. A case of high energy with HAC 
overshoot and spiral shrink can produce both situations. 



For the HAC phase the oscillation was eliminated by reverting back to the acquisition bank 
command and its lower limits when RERRC > RERRLM (7 OOOft) . For the pre-final phase the 
oscillation was eliminated by allowing the bank limit to increase from 30' to 60" as the 
command increases to 100". 

PHILIMIT = MIDVAL(30,60, ,43(ABS(PHIC)) + 17.) 

The bailout mode takes a snapshot of bank, PHIC_ATS, whenever bailout_bank goes from 
false to true, and outputs PHIC_ATS as the bank command. 
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3,13 OVERVIEW OF GRTLS OPEN LOOP GUIDANCE PHASES 6 TO 4 

GRTLS phases 6 and 4 use the angle-of-attack command mode of the flight control system(FCS) 
that is not otherwise used in TAEM. The regular normal acceleration command 
NZC/TAEM/FCS mode is used in phase 5. 

The open loop feature of these phases is that guidance does not close the loop with longitudinal 
state commands to drive ac or NZc, although a and NZ are of course closed loop within the FCS. 
An exception is that phase 4 lateral control is a guidance closed loop direction command that 
is either toward the HAC, or away from it for a high energy S-Tum. 

The objective of this part of GRTLS is to aerodynamical^ reduce altitude rate after falling into 
the atmosphere, and to transition to front side L/D flight at mach 3.2 for normal TAEM closed 
loop guidance energy management. Typical GRTLS atmosphere re-entry profiles show the 
angle-of-attack going from about 10' at ET/SEP to, and holding. 50' for the phase 6 Alpha 
Recovery. 



Mach 

Alpha Profile 


After the rate-of-descent peaks, a calculation is made for the NZ level to transition to phase 5 
NZ Hold at which point a constant NZ near 2 is held to reduce the rate-of-descent. 
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At a rate-of-descent of -HDTRN, 



Rate of Descent Profile 

the transition to phase 4 Alpha Transition, then commands alpha to a profile GRALPR, until 
TAEM phase 1 (or 0) starts at mach 3.2. 

a 14 GRTLS TRANSITIONS - GRTRN 

The transition to phase 5 occurs when NZ builds to a computed value, NZ + DGRNZ. 



Transition Phase 6 to 5 


The alpha profile to be flown in phase 4 is computed in GRTRN.so that a test with it can also be 
used in the transition to phase 4. The transition is mainly on hdot > HDTRN, but also alpha 
must be > GRALPR 
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Mach 


Phase 4 Alpha Transition Profile 


2 


This a Is Increased for phase 4 -Turns, up to the limit of AMAXLD. 

GRALPR = Min(GRALPR / |Cos <{>| , AMAXLD) 


A correction to NZ, 


CORN? NZ - Cos9 / Cos<t>) ( TAS/(T AS+VCO) 

Involves subtracting the C0/C<|) t 'rom the total NZ to get an initial NZC reset of the filter 
used in TGNZC. This replaces the el of NZC last used In phase 5 which would produce a big 
transient at phase 1 initiate if left 

Tie logic for transition to phase IS ." -i (1 or 0) involves energy, NZ, and mach, but the I-Load 
usage of EOWLl=0,and MSW2=MSW1=3.2 produces a transition based only on mach of 3.2. 

Low energy alerts are issued during phase 4 to the pilot for his action. An energy lower than 
EMOH during overhead approach is suggesting to the pilot that he consider downmoding to 
straight in. 



An energy lower than EMEP is advising that the HAC be moved closer to the runway. 
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a 15 GRTLS BODY VERTICAL ACCELERATION - GRNZC 

The NZC for phase 5 is generated as a function of a linear and an exponential term to go from 
the initial to the final constant command value. 




Phase 5 NZC 


a 16 GRTLS ALPHA RECOVERY - GRALPC 

The angle- of- attack command for phase 6 is ALPREC 50". 

The maximum descent rate is captured here, and the DGRNZ used in GRTRN & GRNZC are 
calculated. 
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hdot 


DGRNZ Calculation 

DGRNZ=MIDVAL((HDNOM-HDMAX)DHDNZ,DHDLL,DHDUL) 

-1558 .002 2 -.3 


The angle-of-attack command for phase 4 Is initialized at the a at the start of phase 4. The 
command is then ramped down toward GRALPR at a limit change of GRAL, ± 1 . When a falls 
below GRALPR, then ac Is clamped to GRALPR 



Mach 

Phase 4 Alpha 


3.17 GRTLS SPEEDBRAKE - GRSBC 

The speedbrakes stay tucked in at zero until pitch Jets are shut off at qbar of 20. 


QBAH 



GRTLS Speedbrake 

The command then ramps up to 80.6 and, then at mach 4, ramps back down to 65. for TAEM 
interface. 
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3.18 GRTLS BANK - GRPHIC 


The bank command Is zero until phase 4 at which point the command Is the same as phase 1 In 
TGPHIC to acquire the HAC, and the same as phase 0 In TGTRAN for high energy S-Tums. The 
flag ISTP4 tracks the S-Tum status so that the transition In GRTRN to TAEM will be either to 
S-Tum phase 0, or acquisition phase 1 . 



ES 



S-Tum Energy Profiles 

GRTLS Phase 4 BankCommands 

The only exceptions from normal TAEM operation are that mach must be less than MSW3(7.0), 
and the I-Load PSSTRN is set high at 1000. to enable S-Tums for all overhead approaches. 
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